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Deeply inelastic scattering of electrons oflt nuclei can determine whether parton distributions 
saturate at HERA energies. If so, this phenomenon will also tell us a great deal about how 
particles are produced, and whether they equilibrate, in high energy nuclear collisions. 



1 Introduction 



The deeply inelastic scattering (DIS) experiments at HERA revealed that structure functions grow 
rapidly at small x and large Q^El. At some Xcrit, for a fixed Q^, it is expected that parton distributions 
will saturate-leading to a weaker growth of the structure functions 0. Asymptotically, at least, one 
expects cross-sections to not grow faster than the logarithm squared of the energyu. 

Recently, there have been hints in the DIS data from ZEUS that saturation may be occuring 
already at HERA energies. The derivative of the structure function F2, with respect to log{Q^), is 
plotted as a function of x (each bin in x being averaged over the acceptance for that particular 
x) eI. The data for this quantity, which at high and relatively large x may be simply related to 
the gluon structure function, show the expected rise for decreasing x (and decreasing mean {Q"^)), but 
then flattens and decreases for very small x (and small (Q^)). An interesting feature of the data is 
that the turnover region lies between (Q^) « 1-10 GeV^. The GRV94 DGLAP fit overshoots the data 
in this region. However, other, more recent parton distribution function sets do fit the data. 

One must warn that the data should be approached with some caution since the averaging proce- 
dure is performed over a wide range of for each bin. However, ZEUS has collected more data in 
the low region which will enable them to perform the average over a smaller range of a. We 
will assume here that atleast the qualitative features of the Caldwell plot will remain unchanged. 

A possible interpretation of the data-often called the "Caldwell Plot" -is that it demonstrates the 
onset of parton saturation i. The caveat about the data aside, if indeed we are seeing saturation at 
small X at HERA, this is very relevant for nuclear collisions at LHC, and to a lesser extent for RHIC. 
A sure-fire way to confirm saturation is to collide electrons with nuclei! i. 

The saturation scale determines the typical intrinsic momenta associated with quanta in the nuclear 
wavefunction. These quanta go on shell in a collision, and eventually produce a large multiplicity of 
particles, mostly pions. At very high energies, the saturation scale is the only scale in the problem, 
and one can estimate that the typical momenta of the particles produced in the collision is at this 
scale. If this momenta is large enough, o ne c an estimate in perturbation theory whether or not they 
equilibrate to form a quark gluon plasma Erliy. 

can be studied 



In the following 



Parton saturation, and its implications for both DIS and heavy ion 
systematically in an effective field theory (EFT) approach to QCD at small xO 
section, we discuss parton saturation in this model, both for quarks and gluons. In section 3, we apply 
this model to nuclear collisions. The very early stages of the nuclear collision, where the non-linearities 
in the fields are large, can be studied in this approach. At late times, the fields linearize. Whether the 
partons, which have "emerged" by these times, equilibrate is an interesting issue, of great relevance 
to the quark gluon plasma community. This is discussed briefly in section 4. 

2 Parton Saturation 

In the inflnite momentum frame 00, the effective action for the soft modes of the gluon fleld 

with longitudinal momenta « (or equivalently x = k^/P^ << 1) can be written in light 
cone gauge = as 

Seff = -jd^x-^Gl,G^'''^'' + j^jd\tdx-p%xt,x-)T,{T''W.^^^[A-]{x-,xt)) 

+ i [ d'^xtdx-F[p''{xt,x~)]. (1) 



Here G^^, is the gluon fleld strength tensor, r" are the SU (A^c) matrices in the adjoint representation, 
and W is the path ordered exponential in the x~^ direction in the adjoint representation of SU{Nc), 



W-oo,oc[^ ]ix ,xt) = Pexp 



ig / dx^A^ [x , xt)r" 



(2) 



The above is a general gauge invariant formlli of the action that was proposed in Ref.lll. 

The effective action in Eq. Q above has a remarkable saddle point solution It is equivalent 
to solving the Yang-Mills equations in the presence of the source J'^'" = 6^~^p'^{xt,x~). One finds a 
solution where = 0, and A' = i^Vd^V^ (for i = 1,2) is a pure gauge field which satisfies the 

equation Di^^ = gp{y,x±). Here Di is the covariant derivative di + VdiV\ y = Vo + ^og{x^ /xq) 
is the space-time rapidity, and yo and xq are the space-time rapidity and the Lorentz contracted 
width, respectively, of the hard partons in the fragmentation reeion. At small x, the space-time and 
momentum space notions of rapidity are used interchangeably t3. The momentum space rapidity is 
defined to be y = yo ~ ln(l/x) where x is Bjorken x and yo is the rapidity in the fragmentation region. 
The solution of the above equation is0 

1 / ig rOdy'^p(y',xt)\ .( ig dy' ^ p{y' ,Xt)\^ 



To compute the classical nuclear gluon distribution function, one needs in general to average over 
the product of the classical fields in Eq. ^ at two space-time points with the weig ht F[p\M. For a large 
nucleus, one may approximate F[p] — > J cPxtdy ^^^-^ Tr (/o^), where /U^ is the color charge squared per 
unit area per unit rapidity. The classical gluon distribution for this Gaussian source is 



dN 1 Cf 



1-^w{-^xIxg(x,\]\\ , (4) 




d'^xt 2'Kas x1 \ I 2aCF 

where Ci;' is the Casimir in the fundamental representation and a is the nuclear cross-section'^. 

For large xt (but smaller that 1/Aqcd, the distribution falls like a power law l/xf-and has a 
I /as dependence! For very small xt, the behavior is the perturbative distribution log{xtAQCD)- The 
scale which determines the cross-over from a logarithmic to a power law distribution is, to follow the 
notation of Muelleri, the saturation scale Qs- Setting xt = 1/Qs and the argument of the exponential 
above to unity, one obtains the relation, 

Q's = ^^xG{x,Q',), (5) 

which, for a particular x, can be solved self-consistently to determine Qs- The value of Qs is approx- 
imately 1 GeV for RHIC energies and 2-3 GeV for LHC energies. To compare with the estimate of 
Gyulassy-McLerranlll, set Qs = Cas^/x^ iii the notation of Ref.0'0 where x = Jy'^ /^^(y) and C ~ 5 
for Qs ~ 1 GeV. 

The model calculation above shows that a) the gluon distribution saturates at some scale Qs » 
AQC£)-the momentum distribution of partons grows only logarithmically for kt < Qs, and b) this 
saturation is seen already at the classical level. The gluon distribution has a tail that goes as 1/A;|, 
so one expects the typical intrinsic momentum of the gluons to be peaked at kt ^ Qs- 

A similar behavior is seen for the structure function Fo. The quark distribution in the classical 
gluon field-and F2, can be computed in the same approach till§, and one obtains 

Q'^aNc sr^ 2 r f^^^oc^ , ( "s^i"^ 2 1 



F,(x,Q2) = ^E4/ ^V/''"^^*^*(l-«-p(-^^?-G(x,^)jj 

X 'KlixtAf) (Az\1 - zfQ" + M|) + Kl{xtAf)A) [z" + (1 - zf)\ , (6) 
where is the electric charge squared of a quark of flavor /, 



" Above, we have re-written the expression for the gluon distribution in Ref.a, using the leading log gluon distribution 
to replace [P' and \o%(xtt^QCD) with the gluon distribution xG{x, -i^) at the scale 1/a;?. 



A'j = Q"^ z{l — z) + M'j , Kq^i are the modified Bessel functions, and Ca is the Casimir in the adjoint 
representation. The above equation is the weh known Glauber expression usually derived Jn the 
rest frame of the nucleus. It is heartening that, with the assumption of Gaussian color chargesli^, the 
formalism of Ref.OJ for structure functions in the infinite momentum frame (which, in principle also 
accomodates non-linear evolution for structure functions) reproduces it. For large Q^, it reduces to 
the standard DGLAP expression, while at small it goes to zero as Q^logiQ"^). One then recovers, 
qualitatively, the shape of the Caldwell plot for dF2/dlog{Q'^). 

One obtains from the above equation, in a manner analogous to Eq. ^ the quark saturation scale 
Q1 by replacing Cp — > Ca in Eq. ^. The relative size of the two saturation scales, glue to quark, is 
therefore determined simply by the ratio of the two Casimirs, Ca/Cf- 

What about quantum corrections to the above quark and gluon distributioris? At the one loop 
level, one gets log(l/x) corrections to the Weizsacker-Williams distribution EUE3'0. However, Mueller 
has argued recently that beyond the one loop level, the distribution has the same form as the as the 
above classical gluon distribution. What does change due to small x evolution is the x dependence of 
the saturation scale u. It will be vervinteresting to see if detailed studies of parton evolution in the 
non-linear region confirm this result E3'EJ. 



3 Parton Production 



In the previous section, we discussed the parton distribution in a single nucleus in an effective field 
theory approach to QCD at small x. Kovner, McLerran, and Weigert applied this approach to study 
nuclear collisions at very high energiesBj. One now has two sources of color charge p{x^), one on each 
light cone, described by the currents J^'" = pj^{xt)5{x^) corresponding to — > oo. The classical 
field describing the small x modes in the EFT is then obtained by solving the Yang-Mills equations 
in the presence of these sources. It can be written as A^^ = 9{x'^)9{x~)A'^ + 0^x^)9{—x~^)A'^ + 
9{—x~)9{x~^)A2. Here A12 have A^ = and the transverse gauge fields A\ 2 are, as discussed in the 
previous section, pure gauges. 

A convenient co-ordinate system to study nuclear collisions is x^ = (r, ri^x_i), where r = ^/2x'^x~ 
and rj = \ log(^). Matching the Yang-Mills equations in the different light cone regions at r = 0, one 
can eliminate singular terms by requiring A^ = A\ + A2 and A'^ = ±x^^[A\,A2] at r = (i = 1,2 
are the two transverse co-ordinates). The sum of two pure gauges in QCD is not a pure gauge. The 
coherence of the Weizsacker- Williams fields is broken, and particles are produced-driven by the above 
commutator term. 

Kovner, McLerran, and Weigert solved the Yang-Mills equations perturbatively to 0{p^). Dia- 
grammatically, this can be represented as qq — > qqg. In Schwinger gauge, = 0, and for kt » asp, 
this is the dominant contribution to gluon production at small x. The result has been computed in 
different ways by several authors Bj;^ and it agrees with the pQCD gluon bremsstrahlung expression 
of Gunion and Bertsch at small xc! . 

To the order 0{p^), the saturation effects discussed in the previous section are not visible. The 
result agrees with the pQCD mini-jet result at small x. They become less reliable below the scale 
kt ~ Qs, where all orders in p become important. Given the distribution in Eq. suggesting that 
most of the partons in the small x component of the nuclear wavefunction have intrinsic momenta 
kt ^ Qs, these higher order effects are important for nuclear collisions. They self-consistently regulate 
the divergence at small kt, leading to infrared safe distributions of produced particles. Saturation 
effects may only be marginal at RHIC, but will be absolutely essential at LHC-where most of the 
produced particles will be semi-hard. [—in 

Recently, Krasnitz and I computed gluon production to all orders in p numerically cIrEa. Gauge 
fields in the forward light cone were obtained by solving Hamilton's equations, for each p configuration, 
on a lattice. Observables were computed by averaging over these configurations with a Gaussian 
weight. Assuming boost invariance, the lattice Hamiltonian is the Kogut-Susskind Hamiltonian in 
2+1-dimensions, coupled to an adjoint scalar field. 



The lattice initial conditions, in analogy with the continuum case, were obtained by matching the 
lattice equations of motion in the different light cone regions at r = 0. Our study was performed on 
two dimensional transverse lattices ranging from 10x10 to 160x160. For simplicity, we studied only 
the SU(2) case. 

For large transverse momenta kt ^ asAi, the results of our simulation agree very well with lattice 
perturbation theory. This is the case for both the momentum and time dependence of the fields. The 
lattice coupling is g^jJ-L. For g^/nL ^ 1, one observes large non-perturbative corrections at smaller 
momenta. At very early times, the field strengths are very large-and a parton model description is 
likely not valid. Since a^/i (or Qs) is the only scale in the problem, one expects fields to linearize on 
the time scale ti ~ l/a^/i. This is seen when we compute er on the lattice. It rises from zero and 
goes to a constant at ^ l/a^/i. At late times t ^ t, this quantity is the energy per unit transverse 
area per unit rapidity of produced particles. It scales as fi^. 

To summarize the above, we now have the technology to systematically compute a range of gauge 
invariant observables in high energy heavy ion collisions. An interesting quantity to compute would be 
the number of particles produced per unit rapidity. Several years ago, Blaizot and Muelleri-assuming 
saturation, computed on fairly general grounds the number of produced gluons per unit rapidity to be 

^ = c2AxG{x,Q^,), (7) 

with c = 1. As Mueller has pointed out recently 0, this coefficient cannot be computed precisely 
without knowing the details of the collision. It is now feasible for us to confirm this formula and 
compute the value of the coefficient c|. We hope to report on it in the near future. 

If saturation does occur at ~ 1 GeV and x ~ 10~^, as the HERA data suggest, then depending 
on the X dependence of the gluon density, saturation may occur already for large nuclei at x ~ 10"^- 
which corresponds to RHIC energies. Modulo what c above is, the multiplicity of produced particles 
at RHIC may tell us something about the saturated gluon density at HERA or vice versa. At LHC, 
the gluon densities would correspond to lower values of x for a fixed than at HERA (hence the 
urgent need for an eA collider!). 

4 Parton Equilibration 

At very early times after the nuclear collision, non-linearities in the Yang-Mills equations are extremely 
important, and the concept of partons is not very meaningful. The non-linearities dissipate-producing 
on shell partons in a time scale r oc 1/asfi- Since this time scale is much less than r = ii, it is likely 
the produced partons will further interact, and perhaps equilibrate. (It is only at this stage that it is 
meaningful to think of a parton cascade.) The further interaction of produced partons is beyond the 
scope of the simulations of Krasnitz and myself. However, the "final" distribution of partons in our 
model can, in principle, be used as the initial conditions for a parton cascade. We should also point the 
reader to the Mueller's recent qualitative study of small angle scattering and the onset of equilibration 
in heavy ion collisions 0. From the theoretical viewpoint, the problem of parton equilibration is an 
extremely difficult one-especially since it is not clear that the problem can be formulated in a gauge 
invariant way. 
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